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Abstract 

We compute the couplings of the zero modes and first excited states of gluons, W's, Z 
gauge bosons, as well as the Higgs, to the zero modes and first excited states of the third 
generation quarks, in an RS Gauge-Higgs unification scenario based on a bulk SO (5) x U(l)x 
gauge symmetry, with gauge and fermion fields propagating in the bulk. Using the parameter 
space consistent with electroweak precision tests and radiative electroweak symmetry breaking, 
we study numerically the dependence of these couplings on the parameters of our model. 
Furthermore, after emphasizing the presence of light excited states of the top quark, which 
couple strongly to the Kaluza Klein gauge bosons, the associated collider phenomenology is 
analyzed. In particular, we concentrate on the possible detection of the first excited state 
of the top, t , which tends to have a higher mass than the ones accessible via regular QCD 
production processes. We stress that the detection of these particles is still possible due to an 
increase in the pair production of t 1 induced by the first excited state of the gluon, G 1 . 



1 Introduction 



Five dimensional (5D) warped extra dimensions provide a very attractive beyond the standard 
model physics scenario, since the Standard Model (SM) weak scale-Planck scale hierarchy may 
be explained in a natural way pQ. The observed light quark and lepton masses, as well as the 
suppression of flavor violating operators is naturally satisfied provided the quark and gauge fields 
propagate in the bulk and the first and second generation quark wave functions are suppressed 
towards the so-called infrared brane (IR brane), where the Higgs is localized and where the natural 
scale of energies is of the order of the weak scale [2] , [3] , [I] . 

The propagation of gauge and fermion fields in the bulk leads to the mixing of zero-modes 
with Kaluza Klein (KK) modes, which induces important tree-level effects on precision electroweak 
observables [5], [6]. This happens specially for gauge bosons and third generation quarks [7]-[TT]. 
The latter tend to be localized close to the IR brane in order to generate the large top-quark mass. 
The introduction of a custodial SU(2) R symmetry together with a discrete left-right symmetry 
leads to reduced corrections to the T parameter [12], p3] and helps protect the bottom-quark 
coupling to the Z gauge boson against large tree-level corrections. The top-bottom doublet may 
then be embedded in a bidoublet of the SU{2) L x SU{2) R group. Still, important corrections to 
the precision electroweak observables subsist at the one loop-level, and agreement with data for 
KK masses at the reach of the LHC may only be obtained in a certain region of fermion mass bulk 
parameters of the third generation quarks [H].[15j. 

The above requirements may be satisfied in a natural way by embedding the Standard Model 
gauge SU{2) L x U{1)y group and the global custodial SU{2) R group into an 5*0(5) x U{l)x gauge 
symmetry group [13]. The 5*0(5) x U(l)x symmetry is broken by boundary conditions at the IR 
brane down to SU(2) L x SU(2) R xU(l)x and to SU{2) l xU{1)y at the ultraviolet brane (UV brane). 
The five dimensional components of the gauge bosons associated with the broken gauge symmetries 
at the IR brane have the proper quantum numbers of the Higgs doublet, leading to a natural 
implementation of the Gauge-Higgs unification mechanism [13] [18]. In Ref. [19] the Coleman- 
Weinberg potential for the Higgs field was studied and its dependence on the five dimensional mass 
parameters was determined. It was shown that the region of parameters consistent with precision 
electroweak observables is in good agreement with that required to obtain the breakdown of the 
electroweak symmetry, with the proper values of the top-quark, bottom-quark and weak gauge 
boson masses. 

One of the most important properties of these type of models is the existence of light excited 
states of the top quark. These quarks are strongly coupled to the gauge bosons KK modes which 
are localized towards the IR brane and are light enough so that the first KK mode of the gluon, G 1 , 
tends to decay into them. This, in turn, leads to a reduced decay branching ratio of G 1 into top- 
quarks. These properties, together with an increase in the width of G 1 make the G 1 detection via 
decay into top quarks more challenging than in the models which have been previously analyzed in 
the literature [20] , [21] , [22] . Moreover, for positive values of the bulk mass parameter associated with 
the multiplet of the left-handed top-bottom doublet which is preferred by flavor constraints |23j, 
the couplings of the left- and right-handed top quarks to G 1 become close to each other, leading to 
a reduced left-right top-quark production asymmetry. A similar effect will be present in the LHC 
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phenomenology of the weak gauge bosons KK modes [23] . 

In this work, we shall concentrate on the production of the first excited state of the top quark 
t , at the LHC We shall first show how to consistently determine the couplings of the quarks to 
gauge bosons in a functional way. These couplings will then be used to compute the decay widths 
and production cross section of the first top quark and gluon KK modes. We shall show that the 
presence of G 1 leads to an important enhancement of the t l production cross section for masses 
beyond the ones that can be tested via direct QCD production [32] [35]. This is very important, 
since agreement with precision electroweak observables tends to be obtained for t 1 masses larger 
than 1 TeV, for which detection of t 1 via QCD production becomes very difficult. We shall also 
show that for large ^-masses the branching ratios for the decay of t 1 into a W + and a bottom- 
quark, a Z and a top-quark, and a top-quark and Higgs boson, are in an approximate 2:1:1 relation, 
as required by the Goldstone equivalence theorem. Moreover, we will show that the technique of 
massive jets [33], [35] becomes very important for the reconstruction of the t 1 modes. 

This article is organized as follows. In section 2 we describe our 5-dimensional model. The 
mass spectrum and other relevant parameters of our model consistent with low energy data were 
investigated in Ref. [19]. In section 3 and 4 we will analyze how to obtain the properly normalized 
wave functions for the gauge fields and the fermions, respectively. In section 5 we derive the various 
couplings necessary to study the collider phenomenology for this model and we numerically study 
the dependence of the calculated couplings on the parameters of our model, as well as the decay 
branching ratios. In Section 6 we discuss the collider phenomenology. We reserve section 7 for our 
conclusions. 



2 5-Dimensional Model 

We are interested in a 5D gauge theory with gauge group SO(5) x U (l)x- The geometry of our space- 
time will be that of RSI [I], with an orbifolded extra spatial dimension in the interval x§ G [0, L\. 
The metric for such a geometry is given by 

ds 2 = a 2 (x^)rj flu dx' 1 dx L ' — dx\ . (1) 

where a(x<~) = e~ kxs . The space spanning the fifth dimension corresponds to a slice of AdSs, with 
branes attached at the two boundary points: x§ = (UV brane) and x$ = L (IR brane). 

We place our gauge fields, Am = A M T a and Bm, in the bulk, where T a are the hermitian 
generators of the fundamental representation of SO(5) and generically Tr[T a .T^] = C(5)5 a,/3 . The 
explicit form of the generators [16] are given in Appendix A. Our fermions ip also live in the bulk, 
and they transform under a representation t a of SO(5). 

The 5D action is 

S 5D = I d 4 x [ L dx^g f--LTr{F MN F MN } - -Lg MN G mn + #r N D N _ M )^ , (2) 

J JO V 4 #5 4 §X / 

where = On — iA%t a — iB^ and g$ and gx are the 5D dimensionful gauge couplings. 



2 



The choice C(5) = 1 is a convenient choice, since it allows us to identify the eigenvalues of our 
generators as the weak isospin, with the four dimensional coupling given by g 2 = g 2 /L. Any other 
choice for C (5) may be absorbed into a redefinition of the gauge fields or the gauge coupling leaving 
the physics unchanged. 

To construct a realistic 4D low energy theory, we will break the 5D SO (5) x U(l)x gauge 
symmetry down to the subgroup 5*0(4) x U(l)y = SU{2) L x SU(2) R x U(i)y on the IR brane and 
to SU(2)l x U(1)y on the UV brane, where Y/2 = T 3r + Qx is the hypercharge and Qx is the 
U(l)x associated charge which is accommodated to obtain the correct hypercharge. We divide the 
generators of SO (5) as follows: the generators of SU(2)l,r are denoted by T aL > R and t aL < R , while 
the generators from the coset 50(5)/50(4) are denoted by T a and t a . 

In order to obtain the correct hypercharge and therefore the right Weinberg angle 9w, we need 
to rotate the fields A 3 £ G SU{2) R and B M e U{l) x PU, 



\ a mJ V s c^J-\B m 



(3) 



c = 95 s = 9x (4) 

y/9s + 9 X ' Vdl + 9x 

The correct Weinberg angle is then given by s| ~ tan 2 9w — (0.23/0.77) ~ 0.2987. We will enforce 
A^ to have even parity, corresponding to the hypercharge gauge boson in the 4D low energy limit. 
From now on we will drop the prime on A' 3n , and it will be understood that or refers to 1r, 2r and 
'3 R . 

To implement the breaking of 50(5) on the two branes as stated above, we impose the following 
boundary conditions on the gauge fields: 



d 5 A a ^' Y = A a f' & = A^' Y = 0, x 5 = (5) 

g 5A a L ,a R ,Y = A *= A WK,Y = Q } X 5 = L . (6) 

As discussed in Ref. [T9], this set of boundary conditions on the gauge fields leads to the iden- 
tification of A\ as the Higgs field with a non- vanishing vacuum expectation value (vev). 

We concentrate on the third generation fermions which are the most important for electroweak 
symmetry breaking (EWSB) and electroweak precision test (EWPT) considerations. The SM 
fermions are embedded in full representations of the bulk gauge group as discussed in [Ti].[l5].[T9]. 
We therefore introduce in the quark sector three 50(5) multiplets per generation as follows: 



1L 



X1L -+ 5/3 2/3 \ ,/ ,N 

t ~ D - |^2r(~'+)5/3 Qr(-,+)2/3) m /, ,\ 

&R ~ 02R - A I , \ n ,d( , \ © M M + > +J2/3 

\X2RK-^+)2/3 Qr{-,+)-1/3J 



(7) 
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Q 



3R 



(xU~,+h/3 &{-+)*/* 

\xt R (-+)2/3 q'k d (- +U/3 



,3R ~ 

®T 1R 



^(-,+)5/3 \ / ^(-,+) 5 /3 

+)-l/3/ \^(+,+)-l/3 y 



where we show the decomposition under S"?/ (2) L x SU(2) R , and explicitly write the £7(1) £Af charges. 
The Qs are bidoublets of SU{2) L x SU(2) R , with SU{2) L acting vertically and SU{2) R acting 
horizontally. T x and T 2 transform as (3,1) and (1,3) under SU{2) L x SU(2) R , respectively, while 
u and u' are SU{2) L x SU{2) R singlets. 

As was done in [19] we introduce mass mixing boundary terms, 

C m = 25(x 5 - L) [u' L M Bl u R + QilM B2 Q 3R + h.c] , (8) 

where Mb 1 and M R2 are dimensionless masses. Since the Higgs mixes, amongst other terms, doublets 
with singlets, we see that with the current parity assignments for the singlet component of the first 
multiplet, we would in principle not have positive parity for this coupling at the IR brane. However, 
as discussed in [19], the switch in the singlet parity to (+, +) is completely equivalent to taking 
M Bi ^l/M Bi . 

We also show the parities on the indicated 4D chirality, where — and + stands for odd and 
even parity conditions and the first and second entries in the bracket correspond to the parities in 
the UV and IR branes respectively. Let us stress that while odd parity is equivalent to a Dirichlet 
boundary condition, the even parity is a linear combination of Neumann and Dirichlet boundary 
conditions, that is determined via the fermion bulk equations of motion. 

The boundary conditions for the opposite chirality fermion multiplet can be read off the ones 
above by a flip in both chirality and boundary condition, (— ,+)l — > (+, — )r for example. In 
the absence of mixing among multiplets satisfying different boundary conditions, the SM fermions 
arise as the zero-modes of the fields obeying (+, +) boundary conditions. The remaining boundary 
conditions are chosen so that SU{2) L x SU{2) R is preserved on the IR brane and so that mass 
mixing terms, necessary to obtain the SM fermion masses after electroweak symmetry breaking, 
can be written on the IR brane. Consistency of the above parity assignments with the original 
orbifold Z 2 symmetry at the IR brane was discussed in Ref. [19J. 



3 Gauge Fields 

Solving the equations of motion in the presence of the Higgs vev. h is complicated, as these mix the 
Neumann and Dirichlet modes. However, 5D gauge symmetry relates these solutions to solutions 
with h = [29J. These solutions, which we generally call Jq{x5,0), are related to the solutions in 
the presence of a Higgs vev., /^(xs, h), via a simple gauge transformations: 

f%(x 5 , h)T a = n-\x 5 , h)fg(x 5 , 0)T a Q(x 5 , h), (9) 
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where h) removes the vev h: 



Q(xs, h) = exp 



-iC h hT 4 



J-r, 



dya 2 (y) 



(10) 



and Ch is the Higgs normalization constant chosen such that the Higgs kinetic term is properly 
normalized, Ch = 95(f cr 2 )~ x l 2 . Following the procedure outlined in [19] leads to the following 
wave functions for the gauge fields when h = 0, consistent with the UV boundary conditions: 



fc,n( x 5,0) = C G , n<aL C(x 5 ,m n ), /c.n^jO) = C G ^aS(x 5 ,m n ) 
f G , n (x 5 ,0) = C G ,n,YC{x 5 ,m n ), f G n n (x 5 ,0) = C G)niaR S(x 5 ,m n ) 



(11) 



where the coefficients Cc,n,a are normalization constants and m n is the particular KK mass under 
consideration. The functions C[x5,m n ] and S[x5,m n ] are the solutions to the equations of motion 
in the case of a vanishing Higgs vev., with the following initial conditions: C(0, z) = 1, C'(0, z) = 0, 
S(0, z) = 0, S"(0, z) = z. They are given by [26], [27], 



7T Z 

C(X 5 ,Z) = ^ a ~ 1 ( x 5) 



Y ° [ k )Jl 



7TZ 



S(x 5 ,z) = ^j-a 1 (x 5 ) 



A[- k )y\ 



(— 



ka(x 5 ) 



- Jo [ j I Yi 



\ka{x 5 )J 

(—)' 

\ka(x 5 )J _ 



(12) 
(13) 



We can now calculate f Gn (x5,h), the wave functions with h ^ 0, using Eq. The explicit 
expressions for f G (xs, h) are given in the appendix. The IR boundary conditions give us a system 
of algebraic equations for the coefficients C Gt n t0l . This system of equations can be broken into four 
subsets of dependent equations according to charge. The directions in internal space 1,1^ and 1r 
and 2,2l and 2^ mix and correspond to the components of W +, ~: 



V2 



cos 



S[L\C G;i + sin *g (S[L}C G , R - C[L]C G , L ) = 



(14) 



2 (C GjiL C'[L] + C GM S'[L}) - V2 sin [^j (e 2kL hC h (C[L}C G , L - S[L\C G>iR ) - 2C G -S'[L}) 



+2 cos 



Xah 

fh 



(e 2kL hS[L]C h C G - + C GH C'\L\ - C GM S'[L\) = 



(15) 



2 (C G , iL C'[L] + C GM S'[L\) + v^sin ^ (e 2kL hC h (C[L}C G , L - S[L]C GM ) - 2C G -S'[L]) 



-2 cos 



Xah 
fh 



(e 2kL hS[L]C h C Gi + C G , iL C'[L) - C GM S'[L]) = 



(16) 



with i — 1,2. 

The neutral gauge bosons, Z and 7 are given by a mixture of 3, 3^, 3^ and Y: 
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V2 



cos 



-y/2 sin 
+2 cos 



Xah 
fh 



Xah 

fh 

Xah 

fh 



S[L]C Gti + sin >f (S[L]c^C g , 3r + C[L] ( S<P C GX - C G , 3 J) = 



—\/2c^ sin 
+2c0 cos 



Xah 
fh 



2 (C G)3i + s^Ccy) C'[L\ + 2c C G , 3fll S'[L] 
(e 2fe ^C, (S[L]c^C g , 3r + C[L] (C g , 3l - S<P C G , Y )) - 2C G ^S'[L\) 

(e 2kL hS[L]C h C G j + (C g , 3l - S<f> C G , Y ) C'[L] - c^C g , 3r S'[L\) = 

2 C<P (C g , 3l - S(P C G , Y ) C'[L] + 2 (1 + 4) C g , 3r S'[L\ 
(e 2kL hC h (S[L] C4> C g , 3r + C[L] (-C g>3l + S<P C GX )) + 2C G ^S'[L]) 

(-e 2kL hS[L]C h C G ^ + (-C g>3l + S(t> C GX ) C'[L] + C<p C g , 3r S'[L]) = 



—\Z2Sff, sin 
+2s^, cos 



Xah 

fh 



Xah 
fh 



2 ( 8 +C GflL + (1 + c 2 ) C GX ) C[L] - 2c <p s <t> C G , 3R S'[L] 

(e 2kL hC h {S[L}c^C Gt3R + C[L] {-C g>3l + S<P C GX )) + 2C G ^S'[L\) 

(-e 2kL hS[L]C h C G ^ + (-C g>3l + s<t> C G , Y ) C'[L] + C4> C g , 3r S'[L\) = 



(17) 



(18) 



(19) 



(20) 



Note that in the above, we have dropped the mass dependence of the gauge boson wave functions, 
but it should be understood that m n is taken to be the appropriate mass for the respective modes. 

The procedure to obtain the normalization coefficients will be to drop one of the equations in 
the above subsets and solve for the normalization coefficients in terms of one of them, which can 
then be found by normalizing the related gauge field wave functions to unity. It is a combination 
of these wave functions which share the same charge that turns out to be the true wave function 
for the particular gauge field under study. For example, explicitly for the W + , dropping the first 
equation, we can solve for C G \ and C Gj1r in terms of C Gj i L : 



C G,l - C G, 



where the "Higgs decay constant" is defined as 



—4 cos 


~X G h~ 

fh 


C'[L]S'[L] + V2C h al 2 h sin 


fh 


(S[L]C'[L] + C[L]S'[L]) 


2S'[L] (c h al 2 hcos 


~X G h~ 

fh 


S[L] + v^sin 


~X G h 

fh 







fh 



1 



dilo dya- 2 (y) 
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(21) 
(22) 

(23) 



and X 2 G = 1/2. 



By normalizing the wave functions which contribute to W + } we find Cg,i l through, 



or equivalently, 



\fk*sM 2 + \fG L (*BM 2 + \fh*(*s,h)\ 2 ) dx 



(24) 



C G ,ifkx 5 ,0)\ 2 + |C Gili /^(x 5 ,0)| 2 + \C G , lR f G R (x 5 ,0)\ 2 ) dx 5 = 1. 



(25) 



It is interesting to note that the explicit Higgs dependence in the normalization equations is always 
canceled. This is true for the gauge bosons as well as the fermions. Therefore the only dependence 
of these on the Higgs vev is through the normalization coefficients. 

The condition that the entire system of equations, Eqs. (|14p -(l20 p . has a solution is only non- 
trivially realized if the determinant is zero. This gives us the spectrum of KK masses and was the 
main focus of Ref. [19J. The normalization coefficients must then be calculated using the proce- 
dure defined above. Once all the normalization constants have been computed, we have all the 
information necessary to calculate the wave functions for the gauge bosons with the appropriate 
masses. However, we would like to point out a subtlety in the neutral sector where after electroweak 
symmetry breaking we obtain the Z gauge boson and the photon. In principle, dropping any of 
the four equations and solving the other three should give us the same result. However, since the 
Higgs does not couple to the photon, the photon does not have any component in the 3 and 3# 
directions. Therefore, dropping either f G (L, h) = or f G R (L,h) = 0, does not lead to the photon 
solution in the massless limit. To compute the coefficients associated with the photon, we need to 
drop the boundary condition for either f G L (L, h) or f G (L, h). In this case the Cc, a are given by: 



Cg,3 - Cg,3 l ■ 



G,3 R 



a 



GY 



C'[L) (clafzsm 


Xh 

fh 


-2C[L)S'[L) cot 


Xh 

Vh_ 


) 


y/2S'[L] (afz cos 


Xh 

fh 


+ C[L]S>[L] 


) 



-a 



G,3 



c^C'[L] 
L S'[L] ' 



a 



G,3 L - 



[a L 2 z cos 


Xh 

fh 


- (l + s 2 , - CSC 2 


Xh 
}fh_ 


) S[L)C'[L) + C[L]S'[L]\ 


s 4> (o, L 2 zcos 


Xh 
fh 


+ C[L)S>[L)) 



(26) 
(27) 
(28) 



where we have used the Wronskian relation, 



S'(x 5 , z)C(x 5 , z) - S(x 5 , z)C'(x 5 , z) = za(x 5 ) 2 . 



(29) 



The consistency of this procedure can be seen by the fact that in this case in the massless limit 
f G (xz,, h), Eq.f lC.4l) . and f G R (x5, h), a combination of Eq.f lC.6l) and Eq.( 1C.7l) . are identically zero as 
we expect for the photon field. 
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4 Fermionic KK profiles 



With the introduction of mass mixing boundary terms, Eq. ([S]), and through the Higgs vev, the 
different multiplets are now related via the equations of motion. The details of the calculations for 
the fermion wave functions and their boundary conditions were discussed thoroughly in Ref. [19|. 
The fermion vector functions in the h = gauge and the IR boundary conditions used to calculate 
the coefficients Cf,% may be determined in a similar way as for the gauge bosons. 

Our three fermion multiplets in the h = gauge can be arranged in vector functions in the 
following way: 



/f,1,l(^5,0) 



Cf,iSmx 
Cf^Sm! 
Cf^S-Mx 
Cf^S-Mx 
Cf^Smx 



Cf,qS-m 2 
Cf$S-M2 

Cf,9S-M 2 

Cf,iqSm 2 



fF,3,R{ X 5,fy 



Cf,xiS- 


-Ms 


Cf,uS- 


-Ms 


Cf,i$S- 


-M 3 


C f 14<S- 


-Ms 


Cf,15^- 


-Ms 


Cf,iqS_ 


-M 3 


Cf,itS- 


-M 3 


Cf,\%S- 


-M 3 


Cf,190- 


-M 3 


Cf,2qS 


M 3 



(30) 



where, as for the gauge bosons, the Cp,i are normalization constants. The opposite chiralities 
have opposite boundary conditions and can be read from these ones. The solutions to the fermion 
equations of motion in the h=0 gauge are written in terms of [26] : 



5 f \ n z 
S m {xb,z) = —a 



x 5 ) 



Ji 



+c 



Yi 



Yi 



+c 



Jl 



+c 



ka(x 5 ) 



(31) 



which has initial conditions Sm(0,z) = 0, d$S±M(0, z) = z. The fermion wave function solutions 
with h — 0, Eq. (l3"0~l) . are related to the functions in Eq. flBTi) via: 



S±m(x 5 , z) 



D ±Mx h 



a 2 (x 5 ) 

„±Mx 5 



S±m(x 5 , z) 



S± m (x 5 ,z) = T 



za(x 5 ) 



d<oS±M{x<o, z). 



(32) 



(33) 



The normalization of the fermion wave functions for the free fermion case contains a non-trivial 
dependence on the warp factor, 



dx 5 a(x 5 ) 3 f l F,m n ,(L,R)( X 5i ®)fF,m m ,(L,R)( X 5, 0) — <5 m ,n. 



(34) 



S 



where here f % Fmn r L mi^s, 0) is proportional to either SM(x5,m n ) or SM(x5,m n ) and in this case, 
m n is the mass spectrum in the free fermion case. The superscript i, which denotes the vector 
component, is not being summed over and we denote the chiral indices as L, R for left-handed 
and right-handed fields. It is therefore convenient to eliminate the dependence on a{x^) of the 
normalization condition by redefining the functions 



S- 



±M 



a(x 5 



-S/2 S _ 



±M, 



S- 



±M 



a(x 5 



-3/2^ 



±M, 



or f j 



F,m n ,(L,R) ~^ «(^5j 



-3/2 f i 

Jr 



F,m n ,(L,R) ■ 



(35) 



In the following, we shall apply this redefinition, which allows for a clearer interpretation of the 
localization of the fermion wave functions in the fifth dimension and maintains the validity of the 
functional relations given in Eq. (1301) . 

We can write the final fermion vector functions in the presence of the vev h using the same 
gauge transformation as for the gauge bosons: 



fF,i,L(x 5 ,h) = Qf F)hL {x 5 ,0) 



(36) 

(37) 
(38) 



Applying the boundary conditions at 25 = L, taking into account the mass mixing terms from 
Eq. (jSJ) we derive the conditions on f F (L, h): 



fk'R + M BJkR = ° 

fkt - MsjWt = 



fF,l,R + M Blf F ,2,R - 
f F ,2,L-M B J FXL = 



fi 

/f,3,L 



F,2,L ' 
5,...,10 





= 



(39) 



where the superscripts denote the vector components. 

Asking that the determinant of this system of equations vanishes so that we get a non-trivial 
solution, we obtain the spectrum of the fermion modes [19]. We note here that the same value 
X F = 1/ v^2 is found in the fermionic for gauge bosons. 

To find the normalization coefficients for the fermion wave functions, we use the same procedure 
as for the gauge bosons: solving the equations for the boundary conditions at X5 = L and then 
normalizing the fermions coupled via the Higgs and boundary mass mixing terms. The equations 
are presented in the Appendix. In the equations as well as in the solutions, the functions are 
evaluated at the IR boundary, X5 = L. For the charge 5/3 fermions, the coefficients are: 
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C F n = C F i — = cos 
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F,15 - Cpi8 



Ms 

i 
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Charge —1/3 fermions: 
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Charge 2/3 fermions: 
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(46) 



(47) 



(48) 



Again, as for the gauge bosons, we can use the normalization conditions to solve for the coeffi- 
cients Cf,i, 



l f?l mn ,(L,R)( X 5, h).f%, mnM {x b , h)dx 5 



(49) 



where the index i takes into account the sum over the three different multiplets and the superscript 
Q denotes the charge of the fermion state. Here m n denotes the masses of the zero and KK modes 
fermions of charge Q. We define the vector fp(x 5 ,h) as the vector /f(^5 ? ^) with all components 
which do not participate in the state mixing set to zero. Using the equations of motion for the 
fermion fields and the boundary conditions associated with them, one can show that once one 
imposes the normalization condition for one chirality, the other one is automatically satisfied. We 
checked the consistency of our analysis by numerically verifying the above statement. 
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5 Couplings 
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Figure 1: gGHtlQsik) vs c i- As c\ grows, the coupling for both chiralities unify. 

Once we have the solutions for all the normalization coefficients, we have the necessary informa- 
tion to compute the various boson-fermion couplings. We will start with the coupling of the first 
excited state of the gluon field G l with the fermions. Since gluons do not interact with the Higgs 
field, their profile on the fifth dimension is given by C[x§, z], where in this case z = mp is the G 1 - 
mass. For G 1 -induced pair production of t lJ s and for backgrounds to this production mechanism, 
we are interested in the g G ifi t i, Qcnt and gent 1 couplings. Since the profiles in the 5th dimension 
of left-handed and right-handed fields differ, this leads to left-handed and right-handed couplings 
which are different. Similar to the description given in Eq.([21) for the electroweak interactions, these 
couplings proceed from the covariant derivatives in the fermion kinetic terms, 

gam = g5 S N G i jf fx^^^t^ 5 '^-^^^ 5 '^^ C[x 5 ,m G i]dx 5 (50) 

where to simplify notation we omit the chiral subindices. Furthermore, in the above expression 
is the 5D strong coupling related to the 4D strong gauge coupling by g s = g 5s /\/L and iV^i is the 
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Figure 2: gcnn 1 1 '9s{k) vs c\. As c\ grows, the left-handed coupling remains constant and the 
right-handed one grows. 



normalization for G 1 , which is given by 




L 



1/2 



G 1 = 



C 2 [x 5 ,m G i}dx 5 



(51) 



To obtain g G ifi t i and g G i{i t we only need to replace m t by the appropriate fermion mass in Eq.(l50l). 

Likewise, we calculate the Higgs-fermion coupling gnu 1 , which will be important when studying 
the decay modes of t 1 to Higgs and tops. We look again at the kinetic term for the fermions. The 
Higgs originates from the zero-mode of the A\ gauge field which enters in the fermion covariant 
derivative. Thus the coupling, which mixes left and right handed fields (left-right or right-left), 
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Figure 3: g Z t l tl ' 9w(k) vs m t \ (TeV), where g w (k) = g^/yL. Notice the l/m t i dependence of g z fi 
takes the form, 

= i — C F ,2,lC* F 2fi J e~ kX5 (ifHfF,5fi{fF,2,l + /F,3,l) _ ^(/f,2,0 + /f,3,o) /h/f,5,1 

+ ifHfF,10fi(fF,7,l + /F,8,l) _ ^(/f,7,0 + /f,8,o)/h/f,10,1 + (/f,16,0 _ /f,19,o) ^^/fT' 1 

, ( f * __ f * \ /g/F,13,l , /j* ,* \ /if /.F, 16,1 /,* ,* \ Mot 

+ WF,16,0 ./>,19,oJ UF,12,0 + 7f,13,oJ 1/^,12,0 + 7f,13,oJ 

(52) 

where /# = e 2fca;5 is the Higgs profile and is the Higgs field normalization. For simplicity we 
have written fF,i,m t — fF,i,o, fF,i,m tl = fF,i,i and we again omit the chiral indices. The subscripts 
and 1 on the C's and the fermion functions denote the fermion mass that should be used, m t or 
m t i respectively. Using relationships between the fermion normalization coefficients given in the 
appendix, this can be written as, 
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gmt 1 = — J dxe kxs (Cf,3,oCf,5,iS-m 1 ,oSm 1 ,i — (Cf,5,oCf,2,i — Cf,2,oCf,5,i) Sm 1 ,oSm 1 ,i 

+ 2 F,7 flC f,W,\S -M 2 fiS m 2 ,i ~ Cf,w,oCfj,iSm 2 ,oS-m 2 ,i^ 

+ i\Pl (c F160 (Cp,i2,i + Cf,i3,i) + Cf,5,oCf,3,i£'mi,0'S'-Mi,i _|_ ((7 F12 

+ Cf,13,o) Cf,16,i) S-M 3 ,oS-M 3 ,1 ) (53) 



In the above, the only explicit dependence of the couplings on the Higgs vev comes through the 
fermion normalization coefficients. 

For the case of gmt or g H fi t i, the Higgs couplings take a particularly simple form, 

9mt — J e kxs (Cf,3,oCf,5,o (^Sm 1 ,o[x5]S-m 1 ,o[ x 5\ — Sm 1 ,o[ x 5}S-m 1 ,o[x5]J 

+ 2C F)7) oC F) io,0 (j>M 2 ,o[ X 5\S-M 2 ,o[ X 5\ ~ Sm 2 A X 5}S-M 2 ,oI X 5\J^ 

= - y £ e2kX5 ( C f,3,oC f , 5 ,o + 2C FX oC F ,io,o) dx 5 (54) 
where we have used the Crowian relation: 

- S M (x 5 ,z)S_ M (x 5 ,z) + S M (x 5 ,z)S- M {x5,z) = l/a(x 5 ) (55) 

In the case of g H fi t i, we obtain a similar result with the subindex replaced by 1. 

Another important coupling which will be of major interest for the reconstruction of the t 1 
invariant mass, as will be seen in the collider section, is the g w+ fi b coupling which again comes from 
the covariant kinetic term of the fermion fields, and is given by: 

9w+fib = ^95 (^2 fp,Tm t i ( X5 ' h )' (/oa^s, h)T a ) .f%* mb (x 5 , h) j dx 5 

i^,2,l ^4,0 / p p^ p p p^ p p p^ p 

2 {JG,8Jf,2,1JF,4,0 - JG,5JF,3,lJFA,0 ~ JG,1J F,5,lJF,4fi 

+ «/g,5/f,12,i/f,17,0 + fG,lfF,13,lfF,U,0 + ^-^/g,5/f,14,0 + ^/g,i/f,17,0^ /f,16,1 
+ «/g,8/f,12,i/f,14,0 + /g,5/f,13,i/f,20,0 + ^-^=/g,5/f,14,0 + ^/g,8/.F,20,0^ /f,19,1^ 

(56) 

where these couplings are left-left and right-right, g 5 is the 5.D weak coupling and for simplicity 
we have written fa a (x5, h) = fc,a- Here a runs over the generators associated with the W boson: 
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1, 1l and 1r. The subindex refers to the bottom quark and the subindex 1 to t 1 . Notice that in 
this case the second multiplet does not contribute to the coupling since it doesn't mix to form the 
bottom 5D profile. The first excited state of the bottom fermion has a mass that puts it out of the 
reach for the LHC 

The g Z tH coupling, which is also important as a decay channel for t , is similarly given by, 

9ziH = 95 j Q (X^S^l ( X5 > (fo( X 5, h ) Ta ) -fl'lmM^ h ) + h - C J dx 5 



CF,2flCp2 1 



X 



3#5 



/g,11 — /g,7 + /g,10 /-F,2,o/f,2,1 + /g,3 1 /f,2,0,/f,5,1 



+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 



{^^j f G ' U f G > 7 ~ /g,1oJ /f,3,o/f,3,1 — /g,3/f,7,o/f,10,1 + /g,3,/f,8,0,/f,10,1 



/g,3./f,2,1 + /g,3/f,3,1 + 



4gx 



3#5 



/g,11 — /g,7 + /g,10 /f,7,o/f,7,1 + 



/g,h/f,5,1 J /p,5,0 + /g,3/f,3,o/f,5,1 
4#X 



-/g,3/f,7,1 + /g,3/f,8,1 + ("j^") /g.1i/f,10,1^ 



3#5 
F,10,0 



/g,11 + /g,7 — /g,10 /f,8,o/f,8 



4<? 



A' 



3fi-5 

4#x 



3fl-5 

4#x 

3#5 
3#5 



)2 2 

/f,12,1 + ~/^/g,3/f,16,1 ~~ ~/^/(?,3/f,19,1 ) /f,12,0 

/G,ll + /G,7 — /G,10 ] /f,13,1 ~~ ^/g,3./f,16,1 + ^/g,3,/f,19,1 ) /p,13,0 

/G,ll/F,16,l + 7^/g,3/f,13,1 _ ^/g,3./f,12,1 J /f,16,0 



V2 



/g,3/f,13,1 /f,19,0 



(57) 



where a runs over 3, 3^, 3#, F and we introduce T y = Qx-f = I- 

The functional dependence of the couplings does not change with KK states. The Standard 
Model couplings may be obtained from the above expressions by replacing the appropriate mass. We 
have evaluated numerically the couplings gatt, 9hu, gwtb and gztt for the would be zero- mode masses 
of all the fields and recovered the Standard Model results, as expected, with slight modifications 
coming from the mixing with the KK states which becomes smaller as h becomes much smaller 
than k, where k = ke~ kL gives the characteristic scale of the KK masses. 



5.1 Couplings in the linear regime 

The linear regime is defined as the parameter space of our model where Xh/fh < 0.35. As discussed 
in Ref. [19], all KK masses become larger than the SM particle masses in this regime. The fermion 
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solution, Eq. (13T]) . for z <C k takes the form: 

S M ~z [ ° a- 1 (x 5 )e- 2Mv dy + 0(z 3 ) (58) 



Keeping only up to linear order in h and z, we can use the last expression in Eq. (154p for g^tt 
and see that in this limiting case this coupling reduces to 

'» ~ T = 75 (59) 

recovering the expression for the top Yukawa coupling, y top . Similarly, we can obtain analytical 
expressions for the g w+1 ri LbL , g z iX L t L an d 9m^ R t h couplings in the low energy limit which will be of 
importance when analyzing branching ratios for the decay of t . Since the mass of t 1 is smaller than 
k, one can use an expansion similar to Eq. (!58|) to obtain the previously mentioned couplings. We 
find the following approximate expressions 

g 2 ( h 




(1 + 2 Cl )(l + 2c 2 )y / fcZ v /(l - 2ci)(3 + 2ci)(3 + 2^) 
+ 2 Cl ) 2 (3 + 2ci) + M 2 Bi (l + 2c 2 ) 2 (3 + 2c 2 )) 



9z» LtL - + 4 ail) 



Notice that the appearance m 9w+ti r i, r 

and g Z ti t of an additional power of g, as well as the 
dependence on h/m t i has to do with mixing of the charge 2/3 states through Yukawa coupling 
interactions proportional to h. These approximate expressions are in good agreement with numerical 
values obtained for the full couplings. 
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The decay widths of G 1 to fermions and the decay widths of t 1 to Wb, Zt and Ht are given by, 



i m G x / A m t (22122 m t 

GHt = 487T V 1 ~~ ( 9GH h t h + fGii fl ^ - (#G^ L + 9 G H R t R ~ §9GH L t L gGH R t R )^T 

V G 1 \ G 



.,2 



(63) 

TnVi = (64) 

,2 ™2™2 ^,4 ^,2 



327T m 2 ^ \ m\ x m\ x 

,2 1 „2 



±f(m t >, mz , m ,)'4 + ( 1 + ^ + g + 4 - 24 - 2^f 

87T \ 2 \ rn^! m^x 

2 m t\ 

6 9ziHt L 9zti R t R ™ Z ^ J (65) 

, , (9 2 Ht i tR + 9 2 Ht i tL ) 2 2 2 

tf* 1 * = o r/W' TO *> 5 — ("V + m t _ "V) + 2 9m R t L 9mU R mm t i 



37rrnfi 



(66) 



where 



1 



/(mi, m2, ^3) = - — Y m i + m 2 + m 3 — C 2,m\m\ — 2m\m\ — 2m|mJ (67) 

We are taking the bottom mass to be negligible with respect to the other masses under consideration 
and we only consider g w +£i LbL since g w +i^ H b H < 9 w +£ L b L - For m & > m t ,m w ,m z ,m H , it follows 
that g z iX R t R ^ 9zt l L t L ^ 9Hfi L t R ^ 9m 1 R t L an d hence the expressions for the decay widths of t 1 
reduce to, 

F g W+t~l L b L m ll / Rft x 

Twtlb = ~^r^ (68) 

3 

T 7fl , = zt * 69 

zt 1 32tt m| V ; 

rati* = QO mti (70) 

OZ7T 

It is easy to see in this limiting case that the Goldstone equivalence theorem applies, by simply 
replacing in Eqs.(l68 l) -( j70l) the analytical expressions for the couplings derived in Eqs. (l6"0~l) - (l6"2]) . 
One obtains a relation between the branching ratios of the decays of t l into W + b, Zt and Ht of 
approximately 2:1:1, implying that the main decay channel for t 1 will be through W + b. 



5.2 Numerical Results 

In figures [H and [2] we have plotted the dependence of go^u and g G ifi t i on the bulk mass parameter 
ci which corresponds to the first multiplet. All the couplings were calculated in the linear regime 
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Figure 4: Branching ratios for the decay of t vs 
m t i (GeV). Notice that the 2:1:1 relations holds 
for large m t i . 



Figure 5: Branching ratios for the decay of G l vs 
m G i (GeV). Notice that G 1 decays mostly to t 1 
pairs. 



for values of the parameter consistent with EWSB and EWPT. In figure [T] we notice that as C\ 
goes to positive values, the left-handed and right-handed couplings of the top quark approach each 
other and take values of the order of two times the strong gauge coupling. Moreover, we can see in 
figure [2] that as C\ becomes positive, g G ifi Rt i grows becoming similar in size to the approximately 
common value of the left and right handed ga^tt- Furthermore, #(jifi £t i remains approximately at a 
constant large value of about 2g G i^ t for positive C\. In figure [3] we plot the behavior of the coupling 
QziH as a function of m t i. As seen in this figure, the major contribution comes from the left-handed 
coupling to the fermions. Furthermore, we notice a dependence of gzt J t L ^ V 771 * 1 ; which can be 

understood by looking at Eq. floTl) and taking into account that m t i oc k. The 9w+{i LbL coupling 
shows a similar behavior. 

Apart from the couplings, as was previously discussed in Refs. [14J,[15J,[I2], a very important 
property of this class of models is that the first excited state of the top quark is light enough 
as to enable the decay of G 1 into t 1 pairs. This property, together with the large coupling of 
G 1 to t^s, implies a large branching ratio of the decay of G 1 into t 1 pairs. In figure H] we show 
the branching ratio for the decay of t l into W + b, Zt and Ht. In this figure we can appreciate the 
appearance of the approximate 2:1:1 relation between the different branching ratios, consistent with 
the Goldstone equivalence theorem at large values for the masses of t 1 . In figure Owe observed, as 
predicted, the first KK excitation of the gluon decays mostly to t 1 pairs. The subsequent reduction 
of the branching ratio of the decay of G 1 into top-quark pairs, as well as large width effects, lead 
to important modifications of the collider phenomenology with respect to what has been studied in 
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the literature [20] , [21] , [22] , as we will discuss in detail in the next section. 

6 Collider Phenomenology 

^From the parameter space allowed in this kind of models we see that the typical values for 
the masses of the first excited gauge bosons and third generation up-type fermion Kaluza-Klein 
states are close to 2.5 k and 0.8-1.1 k, respectively. Values of k > 1.4 TeV lead to consistency 
with precision electroweak observables and radiative electroweak symmetry breaking in the linear 
regime [14j,[15j,[19j. This leads to first excited KK gauge bosons with masses larger than about 
3.5 TeV, whose wave-functions are localized towards the IR brane. Since first and second generation 
fermions in these models are localized towards the UV brane, this implies a reduced coupling of 
these fermions to G 1 , which turns out to be on the order of a fifth of the strong gauge coupling, 
rendering the search for these heavy G 1 states at the LHC challenging due to the relatively low 
production cross section. 

Most of the previous studies of G l searches at the LHC have considered the case in which it 
predominantly decays into top quark pairs [20] , [21] , [22] . As we stressed above, in the class of models 
analyzed in this article, which lead to consistency with precision electroweak data and electroweak 
symmetry breaking, the G 1 decays predominantly into t 1 pairs. Moreover, the total decay width 
of G 1 turns out to be in the range of 30 to 40 percent of the G 1 mass, which makes it a broad 
resonance leading to a harder analysis for the reconstruction of the G 1 mass. Therefore, in this 
work we shall analyze the search for the first excited fermion state of the top quark and concentrate 
on the relatively narrow t 1 mass reconstruction. 

Searches for excited states of the top quark, decaying into third generation fermions and either 
Higgs or gauge bosons have been the subject of intensive recent study in the literature. Assuming 
QCD production of these states, it has been shown that masses up to about 1 TeV may be tested 
at the LHC [32] , [33] , [35] . This range of masses, however, falls short of the necessary one to test 
the models under study, since the t 1 -mass tends to be heavier than 1 TeV. In this work, we will 
make use of the increased production cross section of t 1 pairs induced by the presence of G^'s. 
As shown in [15], there is a strong correlation between the masses of the G 1 and the t 1 states, 
777.(31 ~ (2.2 - 3) m^i . For the G 1 -masses allowed by precision measurements, mgi > 3 TeV, the 
niG\-m t \ correlation leads to t l masses above 1 TeV. In figure El we show the different contributions 
to the t l pair production cross-sections for a wide range of t 1 masses beyond those allowed in our 
model. We fix mgi = 4 TeV which corresponds in our framework to m t i ~ 1.3 - 1.8 TeV. It is 
precisely in this region of t 1 masses where the G 1 induced t l pair production cross-section becomes 
larger than the QCD one. We notice in this figure that QCD production and G 1 induced production 
amplitudes interfere constructively. This can be understood, by noticing that, contrary to the QCD 
case, the coupling of the first excited KK gluon to t 1 pairs comes with an opposite sign to the 
coupling of G 1 to first and second generation quarks. In addition, most t l pairs are produced with 
an invariant mass smaller than m G i , leading to a constructive interference between both amplitudes 
due to the behavior of the gluon and G 1 propagators. . 

With this in mind, we are going to show that it is possible to find a signal for t 1 with a mass 
up to about 1.5 TeV when we include the contribution to its production rate of the decay of a G 1 
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state with a mass of about 4 TeV in the high luminosity LHC era. We also consider a second point 
of 5D parameter space where the G 1 -mass is of the order of 3.5 TeV and the t 1 -mass is of order 
1.25 TeV. We concentrate in a region of parameter space with c\ > where possible flavor changing 
operators are suppressed [23]. Moreover, as shown in figures CD and El the suppressed couplings 
9GH L t L ~ 9c 1 t R t R with respect to <7 G i t i t i lead to a large branching ratio for t 1 production induced 
by G lJ s. This can be appreciated in figure [5) 
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Figure 6: Cross section production at the LHC for mgi = 4.0 TeV with couplings 9cn^ L t\ /9s(k) = 
-5.18 and g G x£ R £,l 9s{k) = -2.77. 



6.1 Signal and Background simulation 

We are interested in the study of a mostly SU(2)l singlet, vector-like fermion state t 1 , associated 
with the first KK resonance of charge 2/3. This fermion state, due to gauge invariance, couples with 
regular SM QCD gluons with a strength given by the strong coupling and as stressed before, its 
signatures have been studied in detail in [32],|33j,[35j, showing a potential reach for these particles 
with a mass up to about m t i ~ 1 TeV. 

As explained in the previous section, in the Gauge-Higgs unification models we are considering, 
the G 1 decays mostly to t s, and t 1 prefers to decay mainly in the W + b channel. Thus, one of the 
best options to study this resonance is by means of the lepton+jets final state channel. 

pp^(g + G 1 ) -> tH 1 -> W + bW-b -> rvbbjj, (I = e, /i) (71) 

The jets include quarks of the first two families and the charge conjugated process is also 
considered. The irreducible backgrounds for this signal are coming from the SM by means of ti 



20 



Cl 


c 2 


c 3 


M Bl 


M B2 


h/(V2f h ) 


m G i 


m t i 


9c 1 tt R 


9GHt L 






0.26 


-0.41 


-0.57 


2.2 


0.4 


0.278 


3915.8 


1470.2 


-2.09 


-2.28 


-2.73 


-5.22 


0.24 


-0.41 


-0.58 


2.3 


0.5 


0.318 


3439.6 


1250.5 


-2.12 


-2.50 


-2.67 


-5.20 



Table 1: Points of parameter space chosen for t 1 detection. All masses are given in GeV and the 
couplings are in units of g s (k). 

production and W/Z + jets, and from the exotic it production induced by G 1 . As has been shown 
in [52], [33], the W/Z + jets background can be reduced to negligible levels by requiring two 6-tags 
and lepton+MET. In this work, we impose the above requirements and therefore concentrate on 
the analysis of the main background from it production coming from both QCD and G 1 . 

We simulate the events for the background and signal at partonic level using the Madgraph- 
Madevent package v4.1.33 [36J, which includes the HELAS [37] subroutine to compute the matrix 
elements. The renormalization scales for the QCD background are set at the top mass, and for the 
backgrounds G 1 — ► it and signal G 1 — ► tH 1 are set at the t 1 -mass. Moreover, we choose the set 
CTEQ51 [38] for the structure functions. In order to reduce the amount of events that are produced 
in the background simulation due to the large QCD cross section, we put a mild cut at the generator 
level for the p t of the lepton. This cut implies a weaker restriction in comparison to the cuts that we 
will implement in the analysis and therefore we do not introduce a bias over the background sample 
by implementing it at the generator level. For the signal, we put no cuts at the generator level and 
we produce almost ten times more events than it is expected from the cross section. Therefore to 
obtain a realistic value in the counting of events, it is necessary to normalize to the correct value 
cr£, where a is the corresponding cross section and £ is the LHC luminosity. This procedure allows 
us to obtain smooth profiles and simulate the results at higher luminosities. 

The partonic level events produced by Madgraph are passed through the Pythia [39J and 
PGS 4 [40] packages included in Madgraph. The Pythia package performs the hadronization process 
and includes initial and final state radiation, where at this stage we use the default parameters that 
are included in Pythia. The PGS package performs a simulation of the LHC ATLAS detector based 
on a fit of the CDF data. We then modify the default parameters for the 6-tag efficiencies in order 
to have a more realistic simulation for the high luminosity LHC era. To account for this, we replace 
the b-tag/mistag efficiencies with the set (1/2,1/10,1/30) for underlying fe's, c's and untagged jets 
(gluons and light quarks), respectively [33]. Furthermore, we set the jet cone reconstruction algo- 
rithm to AR = a/ Arj 2 + A(f) 2 = 0.6 in order to optimize the reconstruction of the W mass from 
very collimated jets [33], [35]. This choice of Ai? enhances the ratio of signal over background due 
to the kinematic differences between the W plus b jet system arising from these two processes. We 
will continue with the discussion of this subject in the analysis section, because it is found to be a 
key tool in the search for the ^-resonance. 

In tabled] we specify the 5D parameters, associated couplings and masses which we will use in the 
calculation of cross sections and ^-mass reconstruction. We will denote the first set of parameters 
as point 1 and the second set as point 2. In table [2] we show the values for the cross sections and 
number of events expected at the parton level for the different processes. One interesting aspect 
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of this study as can be seen in table [2] is the relatively large production cross section of top-quark 
pairs compaired to t 1 pairs from G 1 -induced production. That this occurs in spite of the fact that 
r G i 4i < rgi t i t i is a result of the relatively large width of the G l particle that leads to a strong 
departure from the narrow width approximation. Furthermore, enhancements of PDF's at low x 
add to the width effects leading to a preference in the production of tops proceeding from G lj s at 
relatively low p t . The invariant mass of most of the top-quark pairs produced from the G 1 decays 
is, therefore, much smaller than the G 1 mass. This can be seen clearly in figures [7| and [HI where 
both distributions have been normalized to a total of 200 events. 




m, f [GeV] 

Figure 7: tt invariant mass reconstruction from 
QCD and G 1 induced production for point 1. Dis- 
tributions normalized to 200 events. 




m tf [GeV] 

Figure 8: ti invariant mass reconstruction from 
QCD and G 1 induced production for point 2. Dis- 
tributions normalized to 200 events. 



At this point we have not yet included the contribution of K-factors to the cross section. These 
K-factors take into account NLO effects in the computation of the matrix elements, and their direct 
effect is to enhance the value of the cross section for all processes. In |33J it was shown that the 
signal to background ratio, for a heavy quark mass of ~ 600 GeV, computed with and without these 
K-factors is very similar. Moreover, in the computation of the K-factors obtained in Ref. [32J, the 
enhancement of the NLO cross section with respect to the LO one increases with the mass of the 
heavy quark. It was obtain in Ref. [32] that K^tH 1 ) ^ 1.5 and K(it) m 1.3 for m t i = 1 TeV. Thus, 
we expect that the inclusion of these K-factors will enhance our amount of events and signal to 
background ratio. Even in the case in which the signal to background ratio remains approximately 
constant with the inclusion of K-factors, there is still an increase in the significance of about yK 
that we will take into account at the end of our analysis. 

6.2 Event Selection 

The first selection of the events is done at the level of hadronized particles. In order to study the 
properties of signal and background with sufficient statistics, we require that the hadronized events 
fulfill the following two criteria: 
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Point 1 


Point 2 


Process 




iV u Events 




after cuts 


*[/&] 


N° Events 


iV u after cuts 


G 1 — > tt 


4.12 


1236 




1 


4.43 


443 





9 -> t 1 * 1 


0.23 


70 




5 


0.687 


69 


5 


g + G 1 ^ ti 


3025 


907527 




6 


3085 


308509 


6 


g + G 1 -> t 1 * 1 


0.88 


266 




21 


2.015 


201 


13 



Table 2: Value of the cross section at parton level before cuts for final state l~ubbjj including the 
charge conjugated process, expected number of events, and number of events that pass the selection 
cuts in Eqs. ( 1721) and (ITS"]) , for point 1 at 300 fb^ 1 and point 2 at 100 respectively. 

• Isolated lepton with p t > 20 GeV and \r]\ < 2.5 plus missing energy with p t > 20 GeV. 

• At least three jets with p t > 20 GeV and \rj\ < 2.5, with exactly two bottom-tags. 

The lepton+MET requirements are useful to reduce the backgrounds from QCD jets, though we 
will put a stronger cut on the lepton p t in the final analysis. The condition of exactly two b-tags is 
useful to reduce the backgrounds W/Z + jets although it leads to a strong reduction of the signal 
statistics. 

^From the numbers in Table El it is clear that the cross section of the SM background is much 
bigger than the signal. Our approach to reduce backgrounds to manageable levels follows the one 
in [33]. We will put cuts over the pt of the harder bottom and over H t , the scalar sum of the pts. 
The distribution of these two variables for the background and signal are shown in figures [9TTT21 It 
is interesting to notice that the p t distribution of t 1 without the presence of G l is similar to the 
same distribution including the G l . From this behavior we conclude that, apart from the increase 
in the total production cross section, the analysis of the t 1 signal in this scenario would not differ 
significantly from the previous analysis in the search of heavy singlet quarks. 

6.3 Analysis 

To begin with, we would like to discuss the two different ways to reconstruct the M^-mass from 
the hadronic jets which are used in the literature to find this kind of resonances [32] , [33] , [31] . The 
first method is based in the direct search of the two jets coming from the W decay. The first jet is 
selected as the harder jet and the second one is the jet that forms an invariant mass with the first jet 
close to the V^-mass. Thus, to optimize this procedure it is necessary to use a cone reconstruction 
algorithm with a AR ^ 0.4 to avoid contamination from other jets. This method works very 
well when the events have W jets which are not too energetic, like in the tt background and in 
the signal when the t l masses are much smaller than 1 TeV [32]. The second method consists of 
reconstructing the P^-mass using the harder jet information. This method is based on the kinematic 
behavior of the ultra-energetic jets coming from W's. These jets, coming from highly boosted W's, 
are so collimated that almost all of the VF-jet energy is deposited in a cone of AR = 0.6 around of 
the E t weighted barycenter [33]. In the scenario of heavy t lJ s, like the ones present in our model, 
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Figure 9: p t of the faster tagged bottom for point Figure 10: Scalar sum of the all p t for point 1. 
1. Distribution normalized to 200 events. Distribution normalized to 200 events. 
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Figure 11: p t of the faster tagged bottom for point Figure 12: Scalar sum of the all p t for point 2. 
2. Distribution normalized to 200 events. Distribution normalized to 200 events. 
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this method increases the signal and decreases the background. This is due to the fact that since 
the bottom produced from the SM energetic top-quarks often contaminates the W-jets, then the 
reconstruction of the invariant V^-mass is less efficient for the background and the opposite effect 
works for the signal, where the bottom and P^-jet are not that much collimated because the heavy 
t 1 is produced almost at rest due to the broad G 1 -width effect. 

Figures [T3] - [T5l show the distribution of the reconstructed invariant P^-mass using both methods; 
the plotted events are required to pass the first selection after Pythia-PGS level and pj. ,max ^> 150 
GeV. To make the comparison more transparent we use a cone algorithm with AR = 0.4 for the 
first method and AR = 0.6 for the second one. From these histograms it is clear that the more 
efficient method to reconstruct the invariant mass of the W coming from t 1 is the one that treats 
the W as a massive jet. Furthermore, we can use this behavior in the invariant mass of the faster 
jet to strongly reduce the background from tt. 

Using the information from the p t and the invariant VT-mass distribution for signal and back- 
ground, we select the next set of cuts to do the reconstruction of the t 1 mass distribution. In the 
case of p^ max and H t we use different cuts for point 1 (left column) and point 2 (right column), 
respectively, 

p b,max ^ GeV) 

H t ^ 1800 GeV, 

p le P ton ^ 15Q GeV _ 

(72) 



250 GeV, 
20 GeV, 

50 GeV, (73) 

where my/bi is the invariant mass of the Wb system and by hi with i = 1,2 we denote the two 
bottom-quarks which come from the decay of either the t 1 or top pair system. As was mentioned 
before, the cut in the lepton p t is imposed in order to reduce the it background, but leads to a 
strong reduction of the signal. The cut in p 3 t ,max strongly reduces the background while leaving the 
signal after the first cut almost intact. The last cut strongly reduces the background as well. This 
can be understood since in the t-t background case, at least one of the bottoms tends to reconstruct 
the top-quark mass with the W gauge boson; something that doesn't happen with the signal. The 
relevance of the other cuts has been discussed above. 

We present in figures [16] and [17] the reconstructed invariant mass m t i vs. the number of re- 
constructed t lJ s per 100 GeV for both signal and background at 300 fb~ x of LHC luminosity for 
point 1 and at 100 fb^ 1 for point 2. In both cases we notice an accumulation of events with an 
invariant mass close to the physical mass of the t 1 particles, leading to the first hint of a heavy 
quark with a relevant branching ratio decay in the W + b channel. It is apparent from both figures 



p b,max ^ 

H t > 1900 GeV, 

p le P ton ^ 20Q GeV; 



The remaining cuts are common to both points, 

j,max \ 

Pt > 

\m w -mL\ c 



\mwbi - m t \ 
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Figure 13: Invariant reconstructed W mass using 
the method of only one jet for point 1. Distribu- 
tion normalized to 200 events. 



Figure 14: Invariant reconstructed W mass using 
the method of only one jet for point 2. Distribu- 
tion normalized to 200 events. 
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Figure 15: Invariant reconstructed W mass using the method of two jets, for point 1 (similar results 
for point 2). Distribution normalized to 200 events. 
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Figure 16: Distribution of the t 1 mass recon- 
structed hadronically for point 1. The numbers 



Figure 17: Distribution of the t 1 mass recon- 
structed hadronically for point 2. The numbers 



of events represent the results after 300 fb 1 of of events represent the results after 100 fb 1 of 



LHC luminosity. 



LHC luminosity. 



that the background has a much flatter distribution in the same mass range after cuts. Following 
Ref. [33] we have considered all W b pairs in each event. This leads to an increase in statistics by 
doubling the number of points per event and furthermore a broadening of the reconstructed mass 
distribution. 

In order to distinguish among the bottom-quarks which we use to reconstruct the t 1 resonance, 
we can impose a further requirement. We will choose the bottom which has the biggest AR with 
respect to the W gauge boson. This requirement is imposed because the bottom quark and the 
W's coming from t 1 decays usually have a large angular separation since £ x 's are produced mostly 
at low p t . This will lead to a preference of real reconstructed ^'s compared to fake ones in the 
signal events. We plot in figures [18] and [19] the reconstructed invariant mass m t i vs. the number of 
reconstructed t^s per 100 GeV for both signal and background at 300 fb" 1 of LHC luminosity for 
point 1 and at 100 fb" 1 for point 2, with this new requirement. 

In order to estimate the statistical significance of the excess of signal over background for points 
1 and 2 at the respective luminosities, we take into account the total number of events of signal (S) 
and background (B) obtained after cuts. Assuming a Poisson distribution of signal and background 
events, the statistical significance would be given by S/ y/S + B. Using the numbers in Table [2] 
for points 1 and 2, we estimate a statistical significance S/y/S + B ~ 4.04 a at a luminosity of 
300 fb" 1 and S/y/S + B ~ 2.98 a at a luminosity of 100 fb" 1 , respectively. As mentioned before, 
the inclusion of K-factors will tend to increase the total significance of the signal. Taking values of 
the K-factors of roughly 1.5 for both signal and background, we obtain a 5 o significance for point 
1 at 300 fb" 1 . For point 2, instead, ^-discovery cannot be achieved even after the inclusion of K 
factors at a luminosity of 100 fb" 1 . However, an increase in luminosity up to about 200 fb" 1 will 
be sufficient to discover t 1 at 5 a significance. 
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Figure 18: Distribution of the t 1 mass recon- 
structed hadronically for point 1. The numbers 
of events represent the results after 300 fb~ l of 
LHC luminosity, with the additional requirement 
of AR in the Wb system. 



Figure 19: Distribution of the t 1 mass recon- 
structed hadronically for point 2. The numbers 
of events represent the results after 100 fb -1 of 
LHC luminosity, with the additional requirement 
of AR in the Wb system. 



An improved analysis based on an exhaustive treatment of physics and non-physics backgrounds, 
with further optimization of the cuts, inclusion of NLO QCD corrections, and a more realistic 
detector simulation should be done to confirm these results. Keeping this in mind, our preliminary 
results indicate that, after including K-factors, 3 a evidence for the presence of these particles may 
be found already at 100 fb' 1 (60 fb~ l ) and that 5 a discovery can be achieved with about 300 fb" 1 
(200 fb' 1 ) for points 1 and 2, respectively. Figures [T8l and [T9l clearly show the importance of the 
contribution of G^'s in t 1 production (and detection), by comparing the number of events produced 
by the process g + G 1 — > tH 1 with the ones produced by the process g — > tH 1 . If we only had direct 
QCD production the signal would be overwhelmed by the background. 

In order to explore the last point further, we have done a simulation considering only QCD 
production of t 1 pairs at the LHC with a luminosity of 300 fb~ l . Assuming the biggest branching 
ratio into the W + b channel for a given t 1 -mass as seen in figure 4, (BRyy+ b « 0.45), using the 
massive jet technique and similar cuts and efficiencies as was done in the previous simulations, we 
found that the maximum t 1 -mass which can be reconstructed leading to a 5 a discovery at the 
LHC is m t i « 1.1 TeV. This corresponds to a t 1 production cross section of about 30 fb. Similarly, 
for point 1 we have shown that by considering G 1 -induced production of t 1 in addition to QCD 
production at a luminosity of 300 fb -1 , the LHC will be sensitive to t 1 masses up to about 1.5 TeV. 
This corresponds to a t 1 production cross section of about 15 fb as presented in figure [61 Though 
the cross section for point 1 is half the value of the pure QCD cross section necessary to discover 
a 1.1 TeV t 1 , the difference can be understood by the increase in the efficiency of the cuts at high 
mass. 

We can use the above information to obtain an estimate of the reach of the LHC at 300 fb~ x 
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Figure 20: Curves of constant cross section for QCD in addition of G l decay, in (mGi,m ( i) plane. 

for arbitrary m t i and m G i masses. In figure [20] we show the curves representing the points in the 
(m G i,m t i) plane with constant cross section of either 30 fb (blue dotted line) or 15 fb (red solid 
line). The lower bounds in the figure correspond to phase suppression due to the impossibility of 
producing £ x 's from G 1 -decays. For masses m G i — > +oo, the blue curve approaches the maximum 
value for ^-masses that can be tested by QCD production alone. All points between the upper and 
lower blue and red lines have a cross section larger than 30 fb and 15 fb, respectively. Since the 
efficiency of the cuts tends to increase with large ^-masses, one expects that the LHC will be able 
to explore the whole region bounded by the lower and upper blue curves with m t i > 1.1 TeV. This, 
however, underestimates the full reach of the LHC since it doesn't take into account the increase in 
cut efficiency at large t 1 masses. The real sensitivity curve for the case of G 1 -induced production 
in addition to QCD production, will lie in between the upper red and blue curves, interpolating 
the points (+oo, l.lTeV) and (4TeV, 1.5TeV) in (m G i,m t i) space. The shaded area in Figure [201 
represents the region of allowed G 1 and t 1 masses in the Gauge-Higgs unification models considered 
in the previous sections. We can see from the figure that the LHC at 300 fb^ 1 will be able to probe 
a large region of parameter space for these type of models which are consistent with electroweak 
symmetry breaking and electroweak precision tests. 

6.4 Discussion and Outlook 

Let us close this section by discussing some important points related to the search for the t 1 
resonance. In the analysis above, we have considered a flat 0.5 b-tagging efficiency, independent of 
the p t of the tagged jet. This may be an optimistic assumption, particularly at the early stages of the 
LHC. The search for the heavy t l, s analyzed in this work, however, will require high LHC luminosity. 
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If the LHC b-tagging efficiency at high p t for two b-tagged jets proves to be lower than the one 
assumed in this work, it may be helpful to perform the search without b-tagging requirements. 
This will increase both signal and background, particularly the W plus jets background. The cut 
on the pt of the harder jet and the massive W-jet, and the requirement of a high-p t lepton and 
large missing energy may be sufficient to reduce the background to manageable levels. The results 
of Ref. [3l] seem promising in this direction. 

In this work, we have only analyzed the search for the t x -quark in its W plus b decay mode. 
However, as discussed above, t 1 decays about half of the time into a top quark and either a Higgs 
or a Z-gauge boson. For the large t 1 masses we considered in this work, the decay products will 
be highly boosted and therefore the reconstruction of these particles will be most efficient by using 
the technique of massive jets. It would be useful to analyze the LHC reach in the search for t^'s 
including these alternative decay modes. 

Single t l production is a complimentary production channel to the one discussed in this work and 
has been analyzed in the context of the little Higgs Model [121 S3] • Similar to the Pi 1 production 
induced by G 1 , the single t 1 production cross-section also becomes larger than the QCD t 1 pair 
production cross-section for sufficiently large values of the t 1 masses. The study of Ref. |l2l H3] 
shows that the t 1 single production extends the sensitivity to t 1 masses at the LHC well above those 
reachable via the QCD process. It would be interesting to analyze single t 1 production in detail in 
our specific model. In addition to likely extending the reach for the t 1 mass, this could also provide 
information about the Wt l b coupling. 

One important property of the Gauge-Higgs unification model analyzed in this article is the 
suppression of the Higgs production cross section in both the iy-fusion and gluon fusion modes [151 
ITTl \19\ HH HH HH] • The reduction of the Higgs boson coupling to gluons could imply the presence 
of extra scalar states that mix with the Higgs and/or of extra colored states. In the class of Gauge- 
Higgs unification models under analysis, the t 1 will provide the dominant contribution to the gluon 
fusion amplitude in addition to the top quark. In this sense, the discovery of the heavy t 1 mode, 
and the study of its properties, may be essential in order to understand Higgs boson production at 
the LHC. 

7 Conclusions 

We have analyzed the possible detection at the LHC of excited states of the top quark in Gauge- 
Higgs Unification models in warped extra dimensions that are consistent with radiative electroweak 
symmetry breaking and precision electroweak measurements. In these models, the t 1 tends to be 
strongly coupled to the first gluon excitation G l and sufficiently light, so that it can be produced 
from decays of G l . This reduces the branching ratio of the decay of G 1 into top quarks and enhances 
the width of G l , making its reconstruction from top decays difficult. On the other hand, consistency 
with precision measurements may only be obtained for masses of t 1 larger than about 1 TeV. 

After presenting a consistent functional method for the computation of the couplings of the 
fermion modes to the gauge modes, we computed the relevant decay widths necessary for the 
collider phenomenology analysis. An immediate consequence of the heaviness of the t 1 is that its 
decay branching ratios into W + b, tZ and tH states are in an approximate 2:1:1 relation which 
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becomes a better approximation with increasing t 1 -mass. On the other hand, for the large t 1 masses 
we considered, m t i > 1 TeV, the SM decay products tend to be boosted, making the reconstruction 
of the weak gauge bosons from its decay into quarks very difficult. Instead, as recently emphasized 
in the literature, it is convenient to apply the technique of massive jets, that becomes essential for 
an efficient reconstruction of the t 1 states in our case. 

Our analysis shows that the increase in the production cross section of t 1 pairs, induced by the 
presence of G lJ s in the class of models under analysis, extends the reach of the t 1 searches to masses 
not accessible via direct QCD production. There is a strong correlation between the G 1 and t 1 
masses, and we obtain that for a value of m G i of about 4 TeV and with a high LHC luminosity of 
about 300 fb' 1 , the first KK excitation of the top quark with a mass of 1.5 TeV may be detected. 
The discovery of a t 1 in this mass range, possible due to an increase in its production cross section, 
may point towards the presence of a heavy gauge boson resonance contributing to the t 1 production 
rate. It may be possible to obtain information about the G 1 resonance from the analysis of t l 
production. Moreover, while t 1 induces a suppression of the Higgs production via gluon fusion, the 
study of the additional decay channel of t 1 into Higgs and a top-quark may lead to novel Higgs 
production signatures at the LHC. In summary, the search of a singlet vector-like heavy quark 
pair produced at the LHC, opens the possibility to explore the parameter space of Gauge-Higgs 
unification models consistent with precision electroweak measurements and electroweak symmetry 
breaking. 
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APPENDIX 



A 50(5) Generators 

The generators of the fundamental representation of 5*0(5), with Tr[T a .T^] = C(5)5 alS , are given 
by: 



rpO,L,R 



rpa 



e abc (S^ - 5)8-) ± (%Sj - 5]5f) 



—i 



(A.l) 



where T a (a = 1,2,3,4) and T a ^ (a L;R = 1,2,3) are the generators of 50(5)750(4) and 50(4) 
respectively. 

Define the matrices A and B as follows, which gives us the generators in the basis that T 3l and 
T 3r are diagonal: 



A = 



i 

V2 



(-i-l \ 

-i 1 

% 1 

-i 1 

\ v^y 



(A.2) 



i 



/< 


1 























o \ 








— i 


1 


























— % 


-1 




















—i 


1 






































-1 


i 
































y/2 























1 


% 



































-1 


i 
































y/2 


\o 




















1 


i 


o / 



(A.3) 



The generators used in the text should always be understood to be related to the ones defined 
above as A.T a .A^ for the fundamental representation and B.T a .B^ for the adjoint representation. 



B Gauge and Fermion Wave Functions with h ^ 

It should be understood in the following that all the previously defined functions, S±M t , S±m { , S 
and C are functions of X5 and z, where z = m n and m n is taken to be the mass of the particle of 
interest. The functions 9 Gx and 6 Fx are defined as: 
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Ogx = 



7 Fx 



f2 ' 
Jhx 

^Fhf h 

j-2 ' 

J hx 



where, 



f 2 = 

J hx 



gUZ 5 dya-\y) 

The gauge boson wave functions in the presence of the Higgs vev are as follows: 



fa,a(x 5 ,h) 



(C.l) 
(C.2) 

(C.3) 



cos9 Gx C Gi + ±sm0 Gx C G>lR ) S[x] - ^ sin 9 Gx C g>1l C[x] 

(cos9 Gx C G2 + ^-sm9 Gx C G ^ S[x\ - ^ sin 9 Gx C g>2l C[x] 
cos9 Gx C G ^S[x] + ±sin9 Gx (c^C G)3jl 5[x] + (s^C GX - C G)3 J C[x\) 

S[x]C G i 



(C.4) 



fa,a L ( x 5,h) = 

- I ((1 + cos0 G J C G)li C[x] + (v^sin^C Gii + (1 - cos0 G J C g , 1r ) S[x\) " 

± ((1 + cos0 G J C G , 2i C[x] + (V2sin^C Gi 2 + (1 - cos0 G .) C g , 2r ) S[x}) 

\ (((1 - cos0 G J C4> C g , 3r + V2sm9 Gx C G ^) S[x] 

+ ((l + cos9 Gx )C G ^ + (l-cos9 Gx )s 4> C GX )C[x]) _ 

(C.5) 
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fa,a R {x 5 ,h) 



- I ((1 - cos0 G J C GAl C[x] - {V2sme Gx C Gti - (1 + cos0 G J C7 G>1 J 
| ((1 - cos0 G J C g , 2l C[x] - (y/2sme Gx C GA - (1 + cos0 G J C G , 2 J 



(((1 + cos0 G J c^C G , 3i? - v^sin^C^) 

+ ((1 - COS0 G J C G:3l + (1 + COS0 G J S C G ,y) C[x\) . 



(C.6) 



fa,x(x 5 ,h) = [ c^C GjY C[x\ - s^C G ^ R S[x\ ] 



The fermion wave functions in the presence of the Higgs vev are given by the following: 
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(C.8) 
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fF,2,R( X 5,h) 



\ ({l + cos0 F J C fj S-m 2 - (1 - cos9 Fx ) C F fiS-M 2 - V% sin 9 Fx C F , w S M ^j 
\ ((1 + cos0fJ C Ft8 S„M 2 - (1 - cos0 F J C fj S-m 2 - \/2s\n9 Fx C F , w S M ^) 

C F $S-M 2 

75 sin 6> Fa; (C F)7 + C F)8 ) 5'_ M2 + cos 9 Fx C F:W S M2 



(C.9) 



fF,3,R( X 5,h) = 

cos9 Fx C Ft nS-M a — 75 sin^ (CV,i5 — C^is) 5-^3 
| ((1 + cos F J C F ,i2 - (1 - cos 6> f J C F ,i3 - « sin 6» Fa; (CV.ie - Cf,iq)) •S'-Mg 
|((1 + cos 9 Fx ) C F) i3 - (1 - cos 9 Fx ) C F; i2 - i sin 6» Fa; (C F; i 6 - C Fi i 9 )) 

cos9 Fx C Ft uS-M 3 — 75 sin 9 Fx (c Ft nS-M 3 — Cf,2oSm3^ 

\ ((1 + cos9 Fx ) C F ,i5 + (1 - cosflpj C F ,i8 - iv^sin^C^n) S-Afg 
| ((1 + cos 9 Fx ) C F ,i6 + (1 - cos F J C F ,i9 - i sin 6» Fa; (C F ,i2 + C F ^)) S- Mz 
\ ((1 + cos6> F J C F ^S- Mz + (1 - cos9 Fx ) C F) 2oS Ma - i\/2s\n9 Fx C FM S- M ^) 
\ (iV2sm9 Fx C Ftll + (1 - cos9 Fx ) C Ftl5 + (1 + cos0 F J C F , 18 ) S. Ma 
\ (i sin 0^ (C F) i2 + C F ,i3) + (1 - cos F J C F ,i6 + (1 + cos F J CV.ig) S-Mg 

_ | ^v / 2sin6» Fa; C , F) i4S'_M3 + (1 - cos0 F J C FA7 S- Ms + (1 + cos0 F J C F ^S M ^j 

(CIO) 

The boundary equations separate by charge as said previously. For the 5/3 charge C^i, CV,n, C.f,i5 
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and Cf,i8 is one coupled set given by, 
2C Ftl S Ml [L} + M B2 (2 cos 



Xph 
fh 



2 cos 



Xph 
fh 



2M B2 C F ^S Ml [L}-(\ 
(iv^sin heh. c FiU + (l + cos 
(-iv^sin ^ C Fill + (l-cos 



Cjrn + iv^2 sin 
Cf,ii + sin 



Xph 



fh 



Xph 



Xph 
fh 

Xph 



Xph 
fh 



(C.11) 
(C.12) 



(C.13) 
(C.14) 

equation involving C F6 is independent. The —1/3 charged fermions are also given by two sets: 
Cf,4, Cf,i4 ? C F ,n, Cf,2 which includes the bottom, 



) C F ,15 + (l - cos 
) Cf,15 + + cos 



(Cf,15 — CF,18)j <S , -M 3 [^] — 0? 

) Cf,i8^ S-m 3 [L] = 0, 

)C F , 18 )^m 3 W =0 



Aj? h 

fh 



Xph 
fh 



the equation 

'F,4, Cf,14) C 

2C F , 4 S_ Ml [L] + M B2 (2 cos 

C fm S-m 3 [L] - M B2 C fa S-mA l }) ~ iV^sin 
) Cf^o'S'msW + ^v^sin 



2 (cos 



Xph 
fh 



C^uS'-MaW +iV2sin [^] (c f , 17 S_m 3 [£] - C F ,2o<W£])) = 0, 

(C.15) 

C Fj 2oSm 3 [L] — CF,nS-M 3 [L]^ — 0, 

(C.16) 



Xph 

fh 



COS 



(l + cos [^]) C f ,2oSm 3 [£] + (-iv^sin [^] C F , 14 + (l - cos [^]) C F>17 ) S_ Ms [£] = 0, 

(C.18) 

and one independently CV,9- 

The set for the charge 2/3 fermion involve C F ,2, C F ,3, Cf,5, C F j, C F;8 , C f ,io, C F ,i2, C F ,i3, C f ,w, C F ,i9, 
which includes the top and is given by, 



C.F,14 + + cos 
Cf,14 + ^1 — COS 



Xph 
fh 



) Cf,17^ •S'-Ma [L] — 0, 
j Cf,17 j S'-Mg [£] 



(C.17) 



-(1- 



cos 



M B2 ((l + cos [^] ) C F) i 2 - (l - cos 
i±])c F , 3 S- Ml [L} + ((. 

! C F i2 - (l + cos 



Xph 
fh 

1 + COS 



Cf,i3 + i sin 



Xph 



Xph 
fh 



j Cf,2 — \/2 sin 



M B2 (^1 - cos 
— ^1 + cos 



Xph 

fh 



Xph 
fh 



Xph 
fh 



Xph 
fh 



)Cf,sS_ Mi [L} + ((l-cos 



Cf,i3 — « sin 
AFh ])c F>2 + >/2siii 



(Cf,16 — Cf,19)^ S-m 3 [L] 

Cf,5 j S'mi [£] =0, 

(C.19) 

(Cf,16 — Cf,i9)^ S~m 3 [L] 

Cf,5 ) S'mi [£] = 0, 

(C.20) 



Xph 
fh 
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sin 



Xph 



fh 



Cf^S-mAL] + Mb 1 (Cfj + Cf,8J S-m 2 [L] + Cf^Sm^ [L] 



+V2 



cos 



Xph 
fh 



C F , 5 S Ml [L] + M Bi C f ,xoSmM) = 



(C.21) 



sm 



Xph 



M Bl Cf, 2 S Mi [L] + C F , 3 5_ Ml [L] + (C F ,7 + C F ,s) S„m 2 [L] 



Xph 
fh 



a/2 sin 
\pl sin 



fh 



Xph 



fh 



— COS 

C F ,ioS M2 [L] - [\1 + cos 
C f ,ioSm 2 [L] + ( ( 1 - cos 



{Mb^f^Smx [L] — Cf,wSm 2 [L\) — 0, 



(C.22) 



Xph 

fh 



Xph 

fh 



Cf7 — ( 1 + cos 



Xph 
fh 



Xph 



cf,8 ) S-M 2 m — 0, 



(C.23) 



fh 



Cf,8 i S-m 2 [L] — 0, 



(C.24) 



M B 2 (^1 + COS 

1 + cos 

Jh 



Xph 

fh 

Cp,12 



Cf,2 — V% 



sm 



Xph 



1 — COS 



Xph 

fh 



C F ,5j S Ml [L] - ( 1 - cos 
CV 13 + i sin 



Xph 



Xph 

fh 



M, 



B 2 



1 — COS 



Xph 



C F ,2 + a/2 si 



— ( ( 1 — COS 



Xph 



Xph 

fh 



fh 

Cri2 - ( 1 + cos 



sm 



Xph 



Xph 

fh 



Cf,sS-Mi[L\ 

(CV,ie — Cf,v}) \ S-m 3 [L] = 0, 

(C.25) 

( '/•■..-, ) * Wl [£] - ( 1 + cos [^1 ) C^S-j* [L] 

(C.26) 



c 



F,i3-i sm 



//I 



z sm 



i sm 



fh 



(0,12 + C Fj i 3 ) + (1 + COS 

(0,12 + O13) - (1 - cos 



Xph 

fh 



Xph 

fh 



Cf.Vo + ( 1 - cos 



Cf,w - ( 1 + cos 



Xph 

fh 



Xph 
fh 



Cf,19J S-m 3 [L] — 0, 

(C.27) 

Cf,i9J S-m 3 [L] = 0, 

(C.28) 
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